Voronoi Tessellations and the Cosmic Web: 
Spatial Patterns and Clustering across the Universe 
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Abstract 

The spatial cosmic matter distribution on scales of a few 
up to more than a hundred Megaparse^ displays a salient 
and pervasive foamlike pattern. Voronoi tessellations are 
a versatile and flexible mathematical model for such we- 
blike spatial patterns. They would be the natural asymp- 
totic result of an evolution in which low-density expanding 
void regions dictate the spatial organization of the Mega- 
parsec Universe, while matter assembles in high-density fil- 
amentary and wall-like interstices between the voids. We 
describe the results of ongoing investigations of a variety 
of aspects of cosmologically relevant spatial distributions 
and statistics within the framework of Voronoi tessellations. 
Particularly enticing is the finding of a profound scaling of 
both clustering strength and clustering extent for the dis- 
tribution of tessellation nodes, suggestive for the clustering 
properties of galaxy clusters. Cellular patterns may be the 
source of an intrinsic "geometrically biased" clustering. 



1. Introduction: the Cosmic Web 

Macroscopic patterns in nature are often due the collec- 
tive action of basic, often even simple, physical processes. 
These may yield a surprising array of complex and gen- 
uinely unique physical manifestations. The macroscopic or- 
ganization into complex spatial patterns is one of the most 
striking. The rich morphology of such systems and patterns 
represents a major source of information on the underlying 
physics. This has made them the subject of a major and 
promising area of inquiry. 

One of the most striking examples of a physical system 



'The main measure of length in astronomy is the parsec. Technically 
a parsec is the distance at which we would see the distance Earth-Sun at 
an angle of 1 arcsec. It is equal to 3.262 lightyears = 3.086 X 10^'^km. 
Cosmological distances are substantially larger, so that a Megaparsec (= 
10® pc) is the regular unit of distance. Usually this goes along with h, the 
cosmic expansion rate (Hubble parameter) H in units of 100 km/s/Mpc 
(h ^ 0.71). 



displaying a salient geometrical morphology, and the largest 
in terms of sheer size, is the Universe as a whole. The 
past few decades have revealed that on scales of a few up 
to more than a hundred Megaparsec, galaxies conglomer- 
ate into intriguing weblike patterns that pervade throughout 
the observable cosmos. Revealed through the painstaking 
efforts of redshift survey campaigns, it has completely re- 
vised our view of the matter distribution on these cosmo- 
logical scales. The spatial distribution of galaxies is far 
from homogeneous. Instead, we recognize a weblike ar- 
rangement. Galaxies aggregate in striking geometric pat- 
terns, outlined by huge filamentary and sheetlike structures, 
the sizes of the most conspicuous ones frequently exceeding 
100/i~^Mpc. Within and around these anisotropic features 
we find a variety of density condensations, ranging from 
modest groups of a few galaxies up to massive compact 
clusters of galaxies. The latter represent the most prominent 
density enhancements in our universe and usually mark the 
dense intersections of filaments and sheets. They stand out 
as the most massive and most recently formed (individual) 
objects in the Universe. Complementing this cosmic inven- 
tory leads to the existence of large voids, enormous regions 
with sizes in the range of 20 — 50ft.^^Mpc that are practi- 
cally devoid of any galaxy, usually roundish in shape and 
occupying the major share of space in the Universe. 

Of utmost significance for our inquiry into the issue of 
cosmic structure formation is the fact that the prominent 
structural components of the galaxy distribution - clusters, 
filaments, walls and voids - are not merely randomly and 
independently scattered features. On the contrary, they 
have arranged themselves in a seemingly highly organized 
and structured fashion, woven into an intriguing weblike 
tapestry that permeates the whole of the explored Uni- 
verse. The weblike spatial arrangement of galaxies and 
mass into elongated filaments, sheetlike walls and dense 
compact clusters, the existence of large near-empty void re- 
gions and the hierarchical nature of this mass distribution 
- marked by substructure over a wide range of scales and 
densities - are three major characteristics of what we have 
come to know as the Cosmic Web. Its appearance is most 




Figure 1. Image of the weblike galaxy dis- 
tribution in the SDSS galaxy redshift survey 
(3rd data release, 374,767 galaxies). 



dramatically illustrated by the most recently produced maps 
of the nearby cosmos. The 2dF - two-degree field - Galaxy 
Redshift Survey (e.g. f4l) mapped the spatial distribution 
of nearly 250,000 galaxies in two narrow sections through 
the local Universe, out to a depth of more than 300/i^^Mpc. 
The SDSS survey (fig.IT] see e.g. 1 12 1), is mapping the loca- 
tion of up to a million galaxies, covering nearly a quarter of 
the sky and reaching out to a depth of beyond 500/i^^Mpc. 

2. Voronoi Models & the Cosmic Web 

In the cosmological context Voronoi Tessellations rep- 
resent the Asymptotic Frame for the ultimate matter distri- 
bution distribution in any cosmic structure formation sce- 
nario. The Voronoi tessellation is the skeleton of the cosmic 
matter distribution, identifying the structural frame around 
which matter will gradually assemble during the emergence 
of cosmic structure. 

Voronoi tessellations are a versatile and flexible mathe- 
matical model for foamlike spatial patterns. They would be 
the natural result of an evolution in which expanding voids 



dictate the spatial organization of the Megaparsec Universe, 
with matter assembling in the high-density filamentary and 
wall-like interstices between the voids. According to re- 
cent work this is indeed what may be expected in standard 
scenarios of cosmic structure formation. Voronoi models 
would represent the asymptotic limit for which the popula- 
tion of voids would correspond to one single void size and 
excess expansion rate. In this paper we seek to sketch the 
background and ramifications of this idea. 

The premise is that some primordial cosmic process gen- 
erated a random (Gaussian) density fluctuation field. The 
troughs, minima, in this field will become the centres of ex- 
panding voids. Matter will flow away until it runs into its 
surroundings and encounters similar material flowing out 
of adjacent voids. The spatial distribution of the density 
troughs in the primordial density field is dependent on the 
specific cosmological structure formation scenario at hand. 

Within the cellular Voronoi skeleton the interior of the 
Voronoi cells corresponds to voids. The Voronoi planes 
are identified with walls of galaxies. The edges delineat- 
ing the rim of each wall are identified with the filaments 
in the galaxy distribution. The most outstanding structural 
elements are the vertices, corresponding to the very dense 
compact nodes within the cosmic web, the rich clusters of 
galaxies. In general, what is denoted as a flattened super- 
cluster will consist of an assembly of various connecting 
walls in the Voronoi foam. The elongated superclusters or 
filaments usually consist of a few connected edges. This 
may be clearly appreciated from Fig.|5] 

2.1. Voronoi Virtues 

Cosmologically, the great virtue of the Voronoi foam is 
that it provides a conceptually simple model for a cellular 
or foamlike distribution of galaxies. By using these geo- 
metrically constructed models one is not restricted by the 
resolution or number of particles. A cellular structure can 
be generated over a part of space beyond the reach of any 
N-body experiment. Even though the model does not and 
cannot addres the galaxy distribution on small scales, it is 
nevertheless a useful prescription for the spatial distribution 
of the walls and filaments themselves. In all, it makes the 
Voronoi model particularly suited for studying the proper- 
ties of galaxy clustering in spatial cellular patterns. Its ease 
and versatility of construction, and its flexibility with re- 
spect to defining cosmological parameters, has made it into 
an ideal tool for statistical studies and tests of structure find- 
ing and identification techniques. 

3. Gravitational Instability 

Comprising features on a typical scale of tens of Mega- 
parsec, the cosmic web offers a direct link to the matter dis- 




Figure 2. A full 3-D tessellation comprising 
1000 Voronoi cells/polyhedra generated by 
1000 Poissonian distributed nuclei. Cour- 
tesy: Jacco Dankers 



tribution in the primordial Universe and contains a wealth 
of information on the cosmic structure formation procesfl 
It thus represents a key to unravelling one of the most press- 
ing enigmas in modern astrophysics, the rise of the wealth 
and variety of structure in the present-day Universe from an 
almost perfectly smooth, virtually featureless, pristine cos- 
mos. 

The generally accepted theoretical framework for the 
formation of structure is that of gravitational instability. 
The formation and moulding of structure is ascribed to the 

^According to our latest insights the Universe is 13.7Gyr old. With the 
COBE and WMAP microwave background telescopes we have observed 
tiny disturbances in the early Universe, when the Universe was around 
379.000 yrs old. These perturbations were no larger than a factor 10~^. 



gravitational growth of tiny initial density- and velocity 
deviations from the global cosmic density and expansion. 
Overdense regions will initially expand slightly less rapid 
than the cosmic background, reach a maximum size, turn 
around and ultimately condense - dependent on the scale of 
the density perturbation - into a recognizable astrophysical 
object. Underdense regions, on the other hand, do have a 
gravity deficit with respect to the surrounding Universe. 

Three fundamental aspects of the ensuing nonlinear 
gravitational clustering process determine the morphology 
of the resulting matter distribution. Theirs? is hierarchical 
clustering: the first objects to form are small compact ob- 
jects which subsequently merge with their surroundings into 
ever larger features. The second fundamental aspect con- 
cerns anisotropic gravitational collapse. Aspherical over- 
densities, on any scale and in any scenario, will contract 
such that they become increasingly anisotropic. At first they 
turn into a flattened 'pancake', rapidly followed by contrac- 
tion into an elongated filament, possibly to finally collapse 
to become a galaxy or a cluster. The third manifest feature 
of the Megaparsec Universe is the marked and dominant 
presence of large roundish underdense regions, the voids 
(see fig.O. 

3.1. Voids 

Inspired by early computer calculations, |5 | pointed out 
that for the understanding of the formation of the large co- 
herent patterns pervading the Universe it is worthwhile to 
direct attention to the complementary evolution of under- 
dense regions. The gravity deficit in their interior makes 
them expand with respect to the background Universe. 
Meanwhile we see the matter assembling in ever dense pla- 
nar and/or filamentary interstices. 

They form in and around density troughs in the primor- 
dial density field. Because of their lower interior grav- 
ity they will expand faster than the rest of the Universe, 
while their internal matter density rapidly decreases as mat- 
ter evacuates their interior (see fig. [3]). They evolve in the 
nearly empty void regions with sharply defined boundaries 
marked by filaments and walls. Their essential role in the 
organization of the cosmic matter distribution got recog- 
nized early after their discovery. Recently, their emergence 
and evolution has been explained within the context of hier- 
archical gravitational scenarios 1 10|. 

As the voids expand they start to take up an increasingly 
major fraction of the cosmic volume. |[5] made the addi- 
tional interesting observation that the outward expansion of 
voids is accompanied by a tendency of voids to assume a 
spherical geometry. While the above properties refer to the 
evolution of isolated voids, we know that in reality they will 
be surrounded by other density structures and that they will 
themselves contain a range of smaller scale features and ob- 




Figure 3. Void evolution: two timesteps in the 
evolution of a void region in a N-body com- 
puter simulation of structure formation (in a 
SCDM model). Courtesy: Erwin Platen. 

jects. Recent work by fTU\ managed to describe the corre- 
sponding "Void Sociology" through an analytical descrip- 
tion within the context of the excursion set formalism Q. 
Fig.[3]contains an illustration of the evolving Cosmic Web, 
marked by the surge and development of large cosmic voids. 
The evolution of voids is dictated by two processes: 

• Voids merge with surrounding voids into a larger void. 

• Voids may get destroyed when embedded within a larger 
collapsing overdensity. 

When working out the implication of the precise interplay 
between these two processes, one finds that voids tend to 
define a population marked by a self-similar peaked size 
distribution centered around a characteristic void size: at 
any one cosmic instant the void population is dominated by 



voids of nearly the same size (fig. [3]). This indeed forms a 
telling confirmation of the observed foamlike distribution of 
galaxies with the characteristic pattern punctuated by voids 
whose size is « 10 — 3Qh^^ Mpc. 

3.2. From Voids to Voronoi 

A bold leap leads to a geometrically interesting situation. 
Taking the voids as the dominant dynamical component of 
the Universe, we may think of the Megaparsec scale struc- 
ture as a close packing of spherically expanding regions. 
The asymptotic configuration for a "peaked" void distribu- 
tion is one of a single characteristic void size. Pursuing this 
situation we would find a cosmic matter distribution which 
would be organized by a population of equally sized and 
equally fast expanding, spherically shaped voids. 

Matter will collect at the interstices between the expand- 
ing voids. In this asymptotic limit of the outflow being the 
same in all voids, these interstices are the bisecting planes 
between two neighbouring expansion centres: the walls 
and filaments would be found precisely at the (bisecting) 
midplanes between expanding voids. For any given set 
of expansion centres, or nuclei, the arrangement of these 
planes defines a unique process for the partitioning of 
space: the resulting skeleton of the matter distribution 
would be nothing else than a Voronoi tessellation |T7', ^ 
(see fig. lU. In other words, the ASYMPTOTIC description 
of the cosmic clustering process leads to a geometrical 
configuration that is one of the main concepts in the field of 
stochastic geometry: 

VORONOI TESSELLATIONS 

A cosmological reaUsation of this process is called a 
Voronoi foam ||6l. 

4. Voronoi Clustering Models 

Voronoi Clustering Models are a class of heuristic mod- 
els for cellular distributions of matter which use the Voronoi 
tessellation as the skeleton of the cosmic matter distribution 

It is the stochastic yet non-Poissonian geometrical dis- 
tribution of the walls, filaments and clusters embedded in 
the cosmic web which generates the large-scale clustering 
properties of matter and the related galaxy populations. It 
is precisely this aspect which is modelled in great detail by 
Voronoi tessellation. The small-scale distribution of galax- 
ies, i.e. the distribution within the various components of 
the cosmic skeleton, will involve the complicated details of 
highly nonlinear small-scale interactions of the gravitating 
matter. Ideally, well-defined and elaborate physical mod- 
els and/or N-body computer simulations would fill in this 
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Figure 4. Schematic illustration of the 
Voronoi kinematic model. Courtesy: Miguel 
Aragon-Calvo. 

aspect. In the Voronoi models described here we comple- 
ment the geometrically fixed configuration of the Voronoi 
tessellations with a heuristic prescription for the placing of 
particles/model galaxies within the tessellation. 

The Voronoi Kinematic Voronoi Model is based upon the 
notion that voids play a key organizational role in the de- 
velopment of structure and make the Universe resemble a 
soapsud of expanding bubbles [5 |. It forms an idealized 
and asymptotic description of the outcome of the cosmic 
structure formation process within gravitational instability 
scenarios with voids forming around a dip in the primordial 
density field. This is translated into a scheme for the dis- 
placement of initially randomly distributed galaxies within 
the Voronoi skeleton. Within a void, the mean distance be- 
tween galaxies increases uniformly in the course of time. 
When a galaxy tries to enter an adjacent cell, the velocity 
component perpendicular to the cell wall disappears. There- 
after, the galaxy continues to move within the wall, until it 
tries to enter the next cell; it then loses its velocity compo- 
nent towards that cell, so that the galaxy continues along a 
filament. Finally, it comes to rest in a node, as soon as it 
tries to enter a fourth neighbouring void. 

4.1. Voronoi Kinematic Model 

All different Voronoi models are based upon the displace- 
ment of a sample of N Voronoi galaxies. The initial spatial 
distribution of these N galaxies within the sample volume 
V is purely random, their initial locations x„o defined by a 
homogeneous Poisson process. A set of M nuclei or expan- 
sion centres within the volume V corresponds to the cell 
centres, or expansion centres driving the evolving matter 



distribution. 

The first stage of the procedure consists of the generation of 
initial conditions of the Voronoi galaxy distribution, along 
with the specification of the properties (width) of the struc- 
tural elements: 

• Distribution of M nuclei, expansion centres, within the 
simulation volume V . The location of nucleus m is y,„. 

• Generate N model galaxies whose initial locations, x„o 
{n = 1 , . . . , N) are randomly distributed throughout the 
sample volume V . 

• Of each model galaxy n determine the Voronoi cell Va in 
which it is located, ie. determine the closest nucleus ja- 

• The width of the walls, filaments and vertices is set by as- 
suming them to have a Gaussian radial density distribution, 
with wall width i?vv> filament width R-p and vertex width 
i?v input parameters of the Voronoi model. 

The second stage of the procedure consists of the calcu- 
lation of the complete Voronoi track for each galaxy n ~ 
1, . . . ,N. Figure |4] contains a sketch of a typical Voronoi 
galaxy track: 

• The first step is the calculation of the galaxy tracks is the 
dtermination for each galaxy n the Voronoi cell Vq in which 
it is initially located, ie. finding the nucleus ja which is 
closest to the galaxies' initial position x„o- 

• Subsequently, the galaxy follows a radial path within 
Voronoi cell Va- It moves along the path emanating from 
its expansion centre ja, ie. along the direction marked by 
the unity vector e„Q, defined by the line starting at and 
pointing radially outward from jet, r„Q, = ya + RnGna- 
This is pursued until the galaxy intersects the Voronoi wall 

• After this the galaxy's displacement is restricted to the 
wall as the displacement component perpendicular to 
the wall is damped. The galaxy moves towards the one edge 
Aai3j defining the nearest intersection with wall path r„Q,^. 
This continues until it intersects Voronoi edge Kap-y- 

• Subsequently, it pursues its path along the spine of the 
edge Kajs-y as the velocities' component perpendicular to 
the edge has been suppressed ("the self-gravity of the fil- 
ament damps its corresponding velocity"). This continues 
till the galaxy finally arrives at vertex 'E^ap^s- The vertex is 
the final destination of the galaxy, the deep potential well 
of the cluster at its location will keep the galaxy within its 
reach. 

By determining the total Voronoi track for each galaxy 
n it is rather straightforward to compute the location of the 
galaxy at any cosmic epoch t by determining the displace- 
ment x„ that each galaxy has traversed along its path. For 
each cosmic epoch we may therefore easily compute the 
corresponding mass distribution. 



Figure 5. A sequel of three consecutive timesteps within the l<inematic Voronoi cell formation pro- 
cess. The depicted boxes have a size of lOO/i~^Mpc. Within these cubic volumes some 64 Voronoi 
cells with a typical size of 25/i^Mpc delineate the cosmic framework with 32000 galaxies. 



4.2. Kinematic Evolution 

The evolutionary progression within the Voronoi kine- 
matic scheme, from an almost featureless random distribu- 
tion, via a wall-like and filamentary morphology towards a 
distribution in which matter ultimately aggregates into con- 
spicuous compact cluster-like clumps can be readily appre- 
ciated from the sequence of 6 cubic 3-D particle distribu- 
tions in Figure|5][14|. 

The steadily increasing contrast of the various structural 
features is accompanied by a gradual shift in topological 
nature of the distribution. The virtually uniform particle 
distribution at the beginning (upper lefthand frame) ulti- 
mately unfolds into the highly clumped distribution in the 
lower righthand frame. At first only a faint imprint of den- 
sity enhancements and depressions can be discerned. In the 
subsequent first stage of nonlinear evolution we see a devel- 
opment of the matter distribution towards a wall-dominated 
foam. The contrast of the walls with respect to the gen- 
eral field population is rather moderate (see e.g. second 
frame), and most obviously discernable by tracing the sites 
where the walls intersect and the galaxy density is slightly 
enhanced. The ensuing frames depict the gradual progres- 
sion via a wall-like through a filamentary towards an ulti- 
mate cluster-dominated matter distribution. By then nearly 
all matter has streamed into the nodal sites of the cellular 
network. The initially almost hesitant rise of the clusters 
quickly turns into a strong and incessant growth towards 
their appearance as dense and compact features which ulti- 
mately stand out as the sole dominating element in the cos- 
mic matter distribution (bottom righthand frame). 

5. Superclustering vs. Vertex Clustering 

Upon quantitatively testing the spatial distribution of 
galaxies various interesting properties of clustering within 



Voronoi foams are revealed. The most salient and outstand- 
ing properties concern the distribution of the Voronoi ver- 
tices. A first inspection of the spatial distribution of Voronoi 
vertices (Fig. |6] bottom righthand frame) immediately re- 
veals that it is not a simple random Poisson distribution. 

The impression of strong clustering, on scales smaller 
than or of the order of the cellsize Xj^ is most evidently 
expressed by the corresponding two-point correlation func- 
tion ^(r) (Fig. lU. Not only one finds a clear positive signal 
out to a distance of at least r « 1/4 Ac but also - surprising 
at the time of its finding on the basis of similar computer 
experiments 1 16| - the correlation function appears to be an 
almost perfect power-law, 

nf, V r / 

7 = 1.95; ro«0.3Ac. (1) 

The solid line in the log-log diagram in Fig. |6] represents 
the power-law with these parameters, 7 « 1.95 and ro ~ 
0.3 Ac. 

Beyond this range, the power-law behaviour breaks 
down and, following a gradual decline, the correlation func- 
tion rapidly falls off to a zero value once distances are of the 
order of (half) the cellsize. Assessing the behaviour of ^(r) 
in a linear-linear plot, we get a better idea of its behaviour 
around the zeropoint "correlation length" ra ~ 0.5Ac (bot- 
tom righthand frame fig. |6]l. Beyond the distribution of 
Voronoi vertices is practically uniform. Its only noteworthy 
behaviour is the gradually declining and alternating quasi- 
periodic ringing between positive and negative values simi- 
lar to that we also recognized in the "galaxy" distribution, a 
vague echo of the cellular patterns which the vertices trace 
out. Finally, beyond r w 2Ac any noticeable correlation 
seems to be absent. 

^in the following we use the cellsize Ac, ie. the intranucleus distance, 
as unit of distance for the generated Voronoi foams. 
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Figure 6. Two-point correlation function of (Poisson) Voronoi vertices. Distance in units of inter- 
nucleus distance. Left: log-log plot of C(r). Right: lin-lin plot of f (r). 



The correlation function of Voronoi vertices is a surpris- 
ingly good and solid match to the observed world. Maps 
of the spatial distribution of clusters of galaxies show that 
clusters themsvelves are not Poissonian distributed. On the 
contrary, they are highly clustered and aggregrate into huge 
supercluster complexes ||8] [I]. Such superclusters repre- 
sent moderate density enhancements on a scale of tens of 
Megaparsec, still co-expanding with Hubble flow, be it at a 
slightly decelerated rate. 

It sheds an alternative view on the power-law cluster- 
ing with power law 7 « 2 found in the cluster distribu- 
tion. Also, the observed cluster clustering length Tq ~ 
20h^^ Mpc can be explained within the context of a cellu- 
lar model, suggesting a cellsize of Ac « 70h^^ Mpc as the 
basic scale of the cosmic foam. On the other hand, this also 
reveals a complication. The suggested cell scale seems to be 
well in excess of the 25h^^ — 35h^^ Mpc size of the voids in 
the galaxy distribution. Also, it does not correspond to the 
clustering of objects in the walls and filaments of the same 
tessellation network when tying it to the observed galaxy- 
galaxy correlation. The solution to this dilemma leads to an 
intriguing finding for the scaling of vertex clustering. 

6. Geometric Scaling 

A prominent characteristic of superclustering in the ob- 
servational world is that clustering of clusters is consider- 
ably more pronounced than that of galaxies. This character- 
istic ties in with a generally observed property of clustering 
to be dependent on the nature of objects. Amongst clusters 
we observe a trend for more luminous/massive clusters to 



be more strongly clustered. While the two-point correla- 
tion function ^cc ('') of clusters is consistent with it being a 
scaled version of the power-law galaxy-galaxy correlation 
function, with almost the same slope 7 k, 1.8, the correla- 
tion amplitude and clustering length tq increase as the char- 
acteristic "mass" of the clusters in the sample gets larger 

6.1. Vertex Selections 

The kinematic Voronoi model allows a geometric mod- 
elling of mass-dependent clustering properties by assigning 
a "mass" to each vertex. Brushing crudely over the details 
of the temporal evolution, we may assign each Voronoi ver- 
tex a "mass" by equating this to the total amount of mat- 
ter which according to the "Voronoi streaming" description 
will ultimately flow towards that vertex. The related nuclei 
are the ones that supply the Voronoi vertex with inflowing 
matter Evidently, vertices surrounding large cells are ex- 
pected to be more massive. It is reasonably straightforward 
if cumbersome to calculate the "vertex mass" Aiy by pure 
geometric means. The details, turn out to be challengingly 
complex, as it concerns the (purely geometric) calculation 
of the volume of a non-convex polyhedron centered on the 
Voronoi vertex. 

In our computer experiments liT4l [Tsl we set up realiza- 
tions of a (Poisson) Voronoi foam comprising 1000 cells 
with an average size of 25/i^^Mpc. From the full ver- 
tex distribution we selected subsets of vertices, each subset 
comprising vertices with a progressively higher lower mass 
limit. 

A telling illustration of the significant stronger cluster- 




Figure 7. The two-point correlation for a variety of vertex subsamples, selected on the basis of "rich- 
ness/mass". The largest subsample (with weakest correlation function) contains all vertices, the 
riches only the 2.5% most massive ones. Left: log-log plot of ^(r). Right: lin-lin plot of ^{r). 



ing for more massive vertices is shown in Fig |6] It shows 
the vertex distribution in a central sUce through the full 3-D 
cubic distribution. While the top lefthand panel shows the 
distribution of the full Voronoi vertex distribution, the dis- 
tribution of the 12.5% richest vertices (top righthand panel) 
and the 2.5% richest vertices (lower lefthand panel) is ob- 
viously significantly different. Not only do we observe a 
marked increase in clustering strength as the sample in- 
cludes more massive vertices, also the spatial extent of the 
clustering patterns appears to grow as a function of vertex 
richness. When correcting for the possibly confusing influ- 
ence of the sampling dilution, by sampling an equal number 
of vertices from each "selected" sample, the effect is even 
more prominent. When lifting the central 1/8*'* region out 
of the 20% vertex subsample in the (top righthand) frame 
and sizing it up to the same scale as the full box, we ob- 
serve the similarity in point process between the resulting 
(bottom righthand) distribution and that of the 2.5% sub- 
sample (bottom lefthand). It is spatial self-similarity in its 
purest form ! 

6.2 Vertex Correlation Scaling 

To quantify the impression given by the distribution of 
the biased vertex selections, we analyzed the two-point cor- 
relation function for each vertex sample. We computed ^(r) 
for samples ranging from the complete sample down to the 



ones merely containing the 2.5% most massive one^ 

The surprising finding is that all subsamples of Voronoi 
vertices do retain a two-point correlation function (Fig. |7]l 
displaying the same qualitative behaviour as the ^vv ('') for 
the full unbiased vertex sample (Fig. |6j. Out to a cer- 
tain range it invariably behaves like a power-law (lefthand 
frame), while beyond that range the correlation functions 
all show the decaying oscillatory behaviour that already has 
been encountered in the case of the full sample. Nonethe- 
less, we can immediately infer significant systematic trends. 

• Clustering Strength. The first observation is 
that the amplitude of the correlation functions increases 
monotonously with rising vertex sample richness. Express- 
ing the amplitude in terms of the "clustering length" Tq we 
find a striking and almost perfect linear relation (Fig. |9] left- 
hand frame). The "fractal" clustering scaling description of 
ifTTl . according to which the clustering scale tq scales with 
the average object distance L{r), turns out to be implicit for 
weblike geometries: 

^ec(r) = f3 {L{r)/ry ; L{R) = n'^/^ . (2) 

This finding is in line with studies based on cluster samples 
of rich clusters or selected on the basis of their X-ray emis- 
sion 19] 121 . They do indeed show a trend of an increasing 
clustering strength as the clusters are richer (« massive). 

*We use the average distance (R) = n{R)~^/^ between the sample 
vertices for characterizing the richness of the sample. This is based on A„ 
increasing monotonously with subsample richness. 




Figure 8. Illustration of the 'self-similar' vertex distribution. The more massive vertex samples are (1) 
more strongly clustered and (2) define clustering patterns on a larger spatial scale. When properly 
scaled a pure self-similar pattern emerges. See text for explanation. 



• Spatial Coherence. A second significant observation, 
from the large scale behaviour inferred from the lin-lin 
plot, is that ^vv extends out to larger and larger distances 
as the sample richness is increasing. The oscillatory be- 
haviour is systematically shifting outward for the richer ver- 
tex samples, reflecting their more extended clustering pat- 
terns. Even though the basic cellular pattern has a charac- 
teristic scale of only Ac, the sample of the 5% richest nodes 
sets up coherent patterns at least 2 to 3 times larger (also cf. 
fig. Pursuing the systematics of this geometric cluster- 
ing, Fig.|9]also plots the coherence scale ra. It scales almost 
perfectly linear with average vertex distance \y ! ! 

These finding are in line with observational evidence 
that £,cc{r) extends out considerably larger scales than 
the galaxy-galaxy correlation ^gg, possibly out to 50 — 
100/i/i~^Mpc, even though less evident than the increased 
strength of clustering. The inescapable conclusion is that 



weblike geometries manage to define coherent structures 
significantly larger than their basic size. They induce clus- 
tering in which richer objects not only cluster more strongly, 
but also over an ever larger spatial range. ! ! ! 

• Self-Similarity. Combining the behaviour of cluster- 
ing scale scale and coherence scale ra a remarkable self- 
similar scaling behaviour is revealed. The ratio of corre- 
lation versus clustering length is virtually constant for all 
vertex samples, r^jra k, 1.86. Weblike geometries may in- 
volve intriguing self-similar patterns of clustering. 

• Clustering Slope. K final interesting detail on the ver- 
tex clustering scaling behaviour is that a slight and interest- 
ing trend in the behaviour of power-law slope. The richer 
samples correspond to a tilting of the slope. Fig. |9] shows a 
gradual change from a slope 7 « 1.95 for the full sample, 
to a robust 7 w 1.8 for the selected samples. Note than this 
suggestive value is found in the majority of observed galaxy 
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Figure 9. Scaling of Voronoi vertex two-point correlation function parameters, as function of average 
vertex separation in (mass) selected subsample, Xv/Xc- Left: clustering length ro (lower, ^(ro) = 1.0); 
coherence length ra (top, £,{ra) = 0.0). Centre: ratio clustering-coherence length, ra/rQ. Right: power- 
law slope 7. 



and cluster samples ! 

7 Summary: Self-Similar Cosmic Geometry 

The revealed systematic trends of vertex clustering have 
uncovered a hidden self-similar clustering of vertices. It 
forms a tantalizing indication for the existence of self- 
similar clustering behaviour in spatial patterns with a cel- 
lular or foamlike morphology. It may hint at an intriguing 
and intimate relationship between the cosmic foamlike ge- 
ometry and various measures of clustering in the Universe. 
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